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Abstract. We apply the path-integral method to study the multiple scattering and collective flow in 
intensity interferometry in high-energy heavy-ion collisions. We show that the Glauber model and 
eikonal approximation in an earlier quantum treatment are special examples of the more general 
path-integral method. The multiple scattering and collective flow lead essentially to an initial 
source at a shifted momentum, with a multiple collision absorption factor that depends on the pion 
absorption cross section and a phase factor that depends on the deviations of the in-medium particle 
momenta from their asymptotic values. 



INTRODUCTION 

The intensity interferometry, as first proposed by Hanbury-Brown-Twiss (HBT) to mea- 
sure the angular diameter of a star [1], has been applied to optical coherence [2], sub- 
atomic physics, and nuclear collisions [3]. In HBT measurements in high-energy heavy- 
ion collisions, the initial source particles undergo multiple scattering and collective flow. 
As a result, it is conventionally assumed that the initial chaotic source evolves into a 
different chaotic source distribution at thermal freeze-out which becomes the distribu- 
tion measured in intensity interferometry. This conventional assumption is recently sub- 
ject to question. Because the intensity interferometry is purely a quantum-mechanical 
phenomenon, the multiple scattering and collective flow must be investigated within 
a quantum-mechanical framework [4, 5, 6]. Applying the Glauber multiple scattering 
theory at high energies and the optical model at lower energies, we find that multiple 
scatterings lead essentially to an initial source distribution with absorption [4, 5, 6]. Us- 
ing the Feynman path integral method, we further find that the collective flow leads to 
a phase factor that depends on the deviations of the in-medium particle momenta from 
their asymptotic values. Subsequent work by Kapusta and Li [7] supports qualitatively 
some of the earlier results of Ref. [4]. Following this suggestion of a quantum treat- 
ment, Cramer and his collaborators [8] later considered the effects of produced pions in 
a static optical model potential. The work of Cramer et al. [8] is however incomplete, as 
the effects of particle collisions and the collective flow have not been taken into account. 

In this manuscript, we apply the path-integral method for the quantum treatment of 
the multiple scattering and collective flow in intensity interferometry. We show that 
the Glauber model and eikonal approximation in an earlier quantum treatment of the 
multiple scattering [4, 5] are special examples of the more general path-integral method. 



PION ENVIRONMENT AND ELASTIC SCATTERING 



Before we present a quantum treatment of intensity interferometry, it is illuminating to 
investigate the pion environment after a chaotic pion source is produced in the phase 
transition of a quark-gluon plasma. 

At the phase transition temperature of T ~180 MeV, the average pion density is ~0.3 
pions/fm 3 , the pion energy is ~0.5 GeV, and the average CM. energy, (y/s nK ), in a K-K 
collision is ~0.7 GeV. Thus the pion source is dense and pions are energetic after their 
production. As {yjs nn ) is only slightly lower than the p mass, a substantial fraction of 
K-K collisions will go through the 7=1, p resonance. The width of the p resonance is 
150 MeV in free space. In a thermalized medium at T = 150 MeV, the width increases 
substantially to ~300 MeV [9], and the p resonance mean lifetime in the medium 
becomes ~ li/(300MeV) or 0.67 fm/c. The orbiting time for the p meson is of order 
2kt and 2r ~ 0.77 fm [10]. As the p meson mean lifetime in the medium (~0.67 fm/c) 
is much shorter than its orbiting time (~2.4 fm/c), it is unlikely for an I = 1, k-k pair 
to complete an orbital revolution before the p meson breaks apart into two pions. The 
scattering of the two pions through the intermediary p meson is essentially an elastic 
scattering with an energy-dependent amplitude. Furthermore, the kk — > KK threshold 
is ~1 GeV, substantially greater than {y/s KK ). Hence, chemical reactions of pions in the 
medium are essentially completed very shortly after the phase transition. From the state 
of chemical freeze-out to thermal freeze-out, the scatterings suffered by the pions are 
predominately elastic. 



PATH INTEGRAL METHOD FOR INTENSITY 
INTERFEROMETRY 

As a pion propagates from the state of chemical freeze-out to thermal freeze-out, the 
elastic kk scattering can be described by short-ranged scalar and vector interactions, 

v^\q — qi), where q is the coordinate of the propagating pion and qi is the coordinate 
of a pion in the medium. The propagating pion is also subject to a collective flow 
which can be described by a long-range density-dependent mean-field scalar and vector 

interactions, V^ v \q), as in similar cases in the dynamics of the nuclear fluid [11]. The 
Lagrangian for the propagating pion is given by 

L(q,q) = Lmf (q , q) + L col (q,q), (1) 
where L mf (q, q) = - [m K + (q)} y'T^ + ^ . yM {q) _ y 0(v) {q) ^ (2) 

L C0 MA) = ~V^(q)V^e + <l-y^(q)-V^(q), (3) 

V^\q)= N l!v^\q-qi). (4) 

From the Lagrangian of the propagating pion, one obtains the pion three-momentum, 

p = dL/dq = y[ mn + V^(q)]q + \^(q), (5) 



the pion Hamiltonian, 

H = p = pq-L = y\m n + V W (q)}+ V°M (q) , (6) 

and the pion mass- shell condition, 

[pO_ V 0(v) {q)] 2_ [p 2_ v (v) (?)] 2_ [mn + V (s) {q)] 2 = Q 5 (7) 

where 7=1/ a/1 — q 2 and y ( s ' v ) (g) = V^'^ (g) + V^'^ (q) . For a pion produced at x with 
momentum K to propagate in the pion medium to the thermal freeze-out point Xf and be 
detected at the detecting point x d with momentum k, the probability amplitude is [12] 

K(kx -> fcc d ) = / ^ , (g) 

where j$>q... is the sum over all paths q from* to x c i, and the action S(kx — > kx d ;q) is 

/• /-x d ,(pathtf) 
S(kx kx d ;q) = L(q,q) dt = - p(q)-dq. (9) 

We can separate S(jcc — > kx d ;q) into different contributions, 

S(Kx^kx d ;q) = -k- (x d -x) + 8 m f(Kx,kx d ;q) + 8 C0 i(Kx,kx d ;q), (10) 

rx d ,(pa.lhq) rx f ,(path q) 

8 mf (Kx, kx d ;q) = - [p m f(q )-k]-dq = - / [Pmf(<? ) - k] • dq , (11) 

(path 9 ) /-x/^path^) 
8 C(A {Kx,kx d ;q) = - I p co[ {q )-dq = - / p co i(<? ) , (12) 

Pmf(?) = (y[m* + V®{q)]+V* v \q), y[m K + V m f(q)]q + V^)) , (13) 

pcoi(?) = (r^foR^fo). rn2(?)q+vS(?)). (i4) 

Because of the additivity of the collision potentials in Eq. (4), the phase shift for multiple 
collision 5 co i is a sum of the phase shifts for individual collisions, similar to the case of 
the Glauber wave function in multiple scattering [4, 13], 

N K -\ 

S co i(Kx,kx d ;q) = £ <> col j(Kx,kx d ;q), (15) 
i 

where S co \j(Kx,kx d ;q) is obtained from Eqs. (12) and (14) with the potential v co \{q — qi) 
in place of the total collision potential Vcoi^)- The propagation amplitude is therefore 

K(kx — > kx d ) = J@tqexp{— ik- (x d — x) + i8 m f(Kx,kx d ;q) + i8 C0 \(Kx,kx d ;q)}. (16) 

If one makes the approximation that the dominant contribution to the path integral comes 
from the trajectory along the classical path q c for mean-field motion (which need not be 
a straight line), then the amplitude is approximately 

K(Kx—>kx d ) m exp{-ik ■ (x d - x) + i8 m f(Kx,kx d ;q c ) + i8 co i(Kx,kx d ;q c )} . (17) 

For the propagation of an energetic pion, the phase shifts along a straight-line trajectory 
are just those considered in Ref. [4]. 



TWO-PION CORRELATIONS 



For a pion with momentum ?q produced at xi to propagate to momenta fc, at the detecting 
point Xdi, the amplitude is 

*¥{K iXi -> fcfXd/) = A(jQX / )^^) J fi:(jQX ( - -> fefXdi), (18) 

where A(jqx,-) is the production amplitude, and 0o(x;) is a random and fluctuating 
production phase for the chaotic source. The probability amplitude for the production 
of two identical pions (?Ci, K2) at (xi,x 2 ) to be detected subsequently as k\ at x d \ and fc 2 
at x d2 is 

l^^JCiXl ^ fciXji)^! (K- 2 X 2 ^ ^2^2) + (-^1 <">*2) j- (19) 

The probability P{k\k 2 ) for the detection of two pions with momenta (ki,k 2 ) is the 
absolute square of the sum of the above amplitudes from all x\ and x 2 source points. 
Because of the random and fluctuating phase 0o(x;) for a chaotic source, the absolute 
square of the sum of the amplitudes becomes the sum of the absolute squares of the 
amplitudes [3]. One obtains P{k\k 2 ) = P{ki)P{k 2 )[\ + R(hk 2 )], where P(ki) is the 
probability of detecting a pion of momentum ki, and 

P(h)P(k 2 )R(kik 2 ) = ^A(/c 1 xi)A(/c 2 X2)A(jcix 2 )A(/C2Xi) 

x r x 2 

xK(k\xi -> kix d i)K(K 2 x 2 -> k 2 x d2 )K* (kix 2 -> k\x d \)K* {k 2 x\ -> k 2 x d2 ). (20) 
We can evaluate the product K(kix\ — > £4x^1 ^"'(k^xi — > k 2 x d2 ). It is equal to 

*(>C1JC1 - ta!)^ (JC 2 X1 - fc 2 X d2 ) = e -^-^l-^)+*2-fe 2 -x 2 ) 

x j S>q J <&<j exp{i8n&(K\xi,k\Xdi\q) -i8^(K 2 xi,k 2 x d2 ;q') 

+ i8 C0 i(Kixi,kix d i;q) - i8* ol (K 2 xi,k 2 x d2 ;q')}. (21) 

The real part of the phase difference <5 co i ( Jqxi , kix d i ; q) — 8* oX ( K 2 x\ , k 2 x d2 ; q') is station- 
ary when q — q' = and is random and fluctuating when q — q' ^ 0. The sum over 
exp{z5 C0 i(K"i^i, k\x d \;q) — 8* ol (K 2 x\,k 2 x d2 ;q')} is approximately zero when q — q'^0. 
Thus, the phase factor Qyjp{i8 C0 \{K\x\,k\x d \;q) — 8* oX {K 2 xi,k 2 x d2 ;q')} operationally be- 
haves approximately as S(q — q')exp{—2.ym S C0 \(Kix\ : k\x d i;q')}. Upon integrating 
over S>q' , we obtain 

K(kixi -> tai)^*(K2Xi -> fc 2 x d2 ) ~ e-* 1 '^-* 1 )^-^-^ 

x y ^exp{j'5 mf (jCixi,A:ix rf i;^) -^(k^xi,^^;?) -2J^m 5 C0 i(jCixi,&ix<n ;<?)}. 

In the sum over paths in the above integral, the dominant contribution comes 
from the trajectory q c that minimizes the action difference of the mean field, 
{i8 m f(Kixi,kix d i;q) -i8^{K 2 x\,k 2 x d2 ;q)}. We therefore have 

K(k\x\ -> hx dl )K*(K 2 xi -> fc 2 x rf2 ) « exp{-?li • (x rfl -xi) + i* 2 • {x d2 -x 2 ) 
+i8 mi (Kixi,kix dl ;q c ) - i8^f(K 2 xi,k 2 x d2 ;q c ) -lym 8 C0 i(Kixi,kix dl ;q c )}. (22) 



The other factor in Eq. (20) can be approximated in a similar way. When all the factors 
are collected, we obtain 

2 



R(k u k 2 ) = 



d 4 xe'( kl ~ fe 2)-*+*'0mf (feifc2^;?c)+^coi(fcifc2^;?c)p eff ^ 1 ^ 2)JC -) 



(23) 



y/f(Kix)f{K 2 x) 

where Peff(*l*l,*) = (24) 

«/,(path 

<i>mf(kik2,x;q c ) = / {[PlmfO?) ~ [Plmi (<?) ~ k 2\ } • dq, (25) 

<$>co\{k\k2-,x;q c ) = 8 C0 i(Kix,kix dl ;q c ) - 8* ol (K 2 x,k2X d2 ;qc) ~ 2iJ^m8 C0 i(Kix,kix dl ;q c ), 

and /( joc) is the momentum distribution of the initial chaotic source at x with momentum 
K that evolves asymptotically to k [4]. The collective flow (mean-field interaction) 
'distorts' the initial momentum k into the final detected momentum k. The above results 
from the path-integral method contain those of Ref. [4] as special cases. 
The multiple scattering phase ty co \(k\k 2 ,x;q c ) can be simplified to be [13] 

/•x/,(path q c ) f 
-2J r m8 C0 i(Kix,k 1 x d i;q c ) = - n n (q ) a abs dq , (26) 

J X 

where n K (q') is the pion density at (#') and cr a b s is the pion absorption cross section. 



CONCLUSIONS AND DISCUSSIONS 



In the environment after the phase transition, chemical reactions are essentially com- 
pleted, and the collisions between pions are predominately elastic. The elastic propa- 
gation of a pion can be studied by the Feynman path-integral method. We describe the 
collective flow of the pion by its dynamical motion in a density-dependent long-range 
mean-field potential, and its scatterings with other pions by short-range it-it interactions. 
We find that HBT correlation measurements lead to an effective source distribution that 
depends on the initial source distribution at a shifted momentum, with a multiple scat- 
tering absorption factor e l ^ oX and a collective flow phase factor e l ^ mi . 

The multiple scattering absorption factor depends on the pion absorption cross sec- 
tion. A substantial fraction of the it-it collisions will go through the 7=1, p resonance. 
The width of the p meson increases substantially in the medium [9] and the p meson 
mean lifetime in the medium is much shorter than its orbiting time. The scattering of two 
pions through the intermediary p meson is essentially an elastic scattering and does not 
represent an absorption process in intensity interferometry. Therefore, if one considers 
the n + -K + correlation, a detected k + can be absorbed in its propagation through a pion 
medium only by interacting with a n~ in the n + n~ — > n°n° reaction. The cross section 
c{k + k~ — > x°7t°) is equal to (%n/9k 2 ) sin 2 [<5(7 = 0) — 8(1 = 2)], where k is the magni- 
tude of the pion momentum in the K-K CM. system and 8(1) is the it-it phase shift for 
the state of total isospin 7 [14]. One finds an absorption cross section c(lt + lt~ — > lt°lt°) 
of ~8 mb and a mean absorption path length ~12 fm in a fhermalized pion medium at 
r=180 MeV. The mean absorption path length increases as the temperature decreases. 
The degree of lt + absorption by the it + it~ — > it°lt° reaction is small. 



The collective flow leads to a net phase shift (j>mf(kik2,x;q c ) in Eq. (25) that depends 
on the deviations of the in-medium particle momenta Pi m f(q) from their asymptotic 
values kf. The contributions of the two terms from p\ m f and p2 m f in Eq. (25) tend 
to cancel and give rise only to a small effect, as indicated by similar hydrodynamical 
calculations where one evaluates m f by following pion trajectories [6]. Because of 
the small absorption due to multiple scattering and the small net phase shift due to 
the collective flow, the effective density measured in HBT measurements is expected 
to depend on a source distribution close to the initial (chemical freeze-out) source 
distribution, at a shifted momentum. We expect that as the initial pion source transverse 
dimension is approximately the spatial dimension of the colliding nuclei, the effective 
distribution measured by HBT should be nearly independent of the collision energy. We 
also expect that in the initial source prior to the collective flow, the transverse size in the 
"out" direction should be approximately the same as that in the "side" direction. Hence, 
Rom/Rside should be approximately close to 1. These expectations are consistent with 
the gross features of HBT transverse radii in high-energy heavy-ion collisions. It will 
therefore be of great interest to carry out numerical calculations of the evolution of the 
produced hadron matter to study the above effects. 
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